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The calculation of finite temperature contributions to the scalar potential in a quantum 
field theory is similar to the calculation of loop corrections at zero temperature. In natu¬ 
ral extensions of the Standard Model where loop corrections to the Higgs potential cancel 
between Standard Model degrees of freedom and their symmetry partners, it is interesting 
to contemplate whether finite temperature corrections also cancel, raising the question of 
whether a broken phase of electroweak symmetry may persist at high temperature. It is well 
known that this does not happen in supersymmetric theories because the thermal contribu¬ 
tions of bosons and fermions do not cancel each other. However, for theories with same spin 
partners, the answer is less obvious. Using the Twin Higgs model as a benchmark, we show 
that although thermal corrections do cancel at the level of quadratic divergences, subleading 
corrections still drive the system to a restored phase. We further argue that our conclusions 
generalize to other well-known extensions of the Standard Model where the Higgs is rendered 
natural by being the pseudo-Nambu-Goldstone mode of an approximate global symmetry. 
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I. INTRODUCTION 

Naturalness of the Standard Model (SM) requires the cancellation of divergent contributions 
to the Higgs mass at the loop level. Most known solutions of the little hierarchy problem involve 
introducing new particles that cancel the divergences caused by their SM partners, where the 
cancellation relies on the existence of a symmetry. In the case of supersymmetry (SUSY) the 
symmetry in question is a spacetime symmetry that relates bosons and fermions whereas models 
that realize the Higgs field as a pseudo-Nambu-Goldstone boson (pNGB) accomplish this with an 
internal symmetry. To be precise, in this paper we will assign a very specific meaning to the word 
“natural”, namely we will label a model as natural if (after cancellations) any existing quadratically 
divergent contributions to the Higgs potential are of the same order as, or negligible to the leading 
logarithmic contributions. 

In the SM, electroweak symmetry is restored at temperatures above 0(100) GeV|IH6]. Exten¬ 
sions of the SM display a similar behavior at finite temperature. In particular, finite temperature 
breaks SUSY, and therefore the diagrams whose quadratic divergences cancel each other at zero 
temperature no longer cancel at finite temperature, generating a thermal mass for the Higgs pro¬ 
portional to T 2 and restoring electroweak symmetry. On the other hand, there is no a priori reason 
why the cancellation of quadratic divergences should not persist at finite temperature for models 
with same-spin partners. It was investigated in ref. [ZH2] whether this may lead to the existence of 
a broken phase of electroweak symmetry at high temperature in Little Higgs models |10Hl4Q In 
the following, we will perform a similar analysis on a theory with same-spin partners. 

A calculation keeping only the quadratically divergent but not subleading contributions to the 
Higgs finite temperature effective potential can be justified when SM partners thermally populate 
the plasma, which happens only close to the cutoff of the effective held theory. However, in 
theories such as the Littlest Higgs pH] where the Higgs is nonlinearly realized as a pNGB, higher 
dimensional terms in the effective potential can become important close to the cutoff due to power- 
law divergent contributions. Then, a one-loop analysis may prove insufficient for calculations at 
energies above the decay constant / of the sigma model. In fact, the effective held theory of the 
Littlest Higgs nonlinear sigma model becomes strongly coupled well below 47r/ due to higher¬ 
dimensional operators being corrected by scalar loops [18| . Thus, the analysis in ref. [7j with the 
Littlest Higgs EFT is untrustworthy for T > / since Matsubara modes have masses of order -kT. 
For this reason, we will choose a benchmark model in this paper which has a weakly coupled linear 

1 Of course, the subject of possible symmetry non-restoration has a long history that significantly predates the Little 
Higgs mechanism, starting with refs. I15lil7| . 
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UV completion, namely the Twin Higgs jl9j . where a one-loop calculation should be reliable. 

In our calculation, we will include subleading corrections in the finite temperature potential, 
which are of a size comparable to the zero-temperature effective potential, and therefore cannot 
be neglected. We find that while we agree with ref. [7] that the thermal corrections of O ( T 2 ) 
do cancel, subleading corrections still restore the symmetry at high temperature in the Twin 
Higgs model. Furthermore, we will argue that our conclusions extend beyond the Twin Higgs 
model, and should remain valid in models where the cancellation of O (T 2 ) corrections to the 
Higgs potential are ensured by an approximate global symmetry of which the Higgs is a pNGB, 
and therefore electroweak symmetry is generically restored at high temperature in models that are 
natural according to our definition. 

This paper is organized as follows: We review the salient features of the Twin Higgs model 
in section [TTJ followed by a review of the general aspects of calculating the finite temperature 


effective potential in section III We then calculate the finite temperature effective potential for 
our benchmark model and we present the results in section IV In section [V] we consider the 
symmetry structure of other well-known natural extensions of the SM where the Higgs is realized 
as a pNGB and we conclude that the lessons learned from the benchmark model are generic. 


II. THE TWIN HIGGS MODEL 

There are several variations keeping with the spirit of the Twin Higgs setup |19H29| and here 
we adopt a minimal version of the model presented in ref. (29] as a benchmark model, and limit 
ourselves to a description of the aspects most relevant to our purposes. The reader is invited to 
consult the original references for any additional details not presented here. 

In very rough terms, Twin Higgs models introduce a second set of degrees of freedom identical 
to the SM. The second set of fermion fields fill out the same gauge representations under the 
new gauge groups as the SM fermions do under the SM gauge groups. The two sectors couple 
to each other through the scalars (Higgs), and in our benchmark model, they are both charged 
under C/(l)y. Furthermore, an approximate Z 2 symmetry relates these two sectors (with sector 
A being identified as the SM). Since the main interest in constructing the Twin Higgs setup is 
to keep contributions to the Higgs potential under control, in many phenomenological studies, 
all fermion fields are neglected, for simplicity, except those that are relevant for cancelling the 
divergent contributions due to the top Yukawa coupling, and this is the approach that we adopt 
as well. 
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For the purposes of this study, we will take the gauge symmetry of the theory to be 

G = [517(3) x SU( 2)] 2 x U(l) Y = [SU( 3) x SU{2)} A x [517(3) x SU{2)} B x U( l) y , (1) 


and the relevant fermionic degrees of freedom in the top sector fill out the representations Qa,b = 
( 3,2 ) AB ,n, B — (3, 1)a,b with hypercharges 1/6 and —2/3, respectively. Under the Z 2 symmetry, 
the gauge and matter fields of the A and B sectors are exchanged (and the U( l)y is unaffected). 

It should be noted that this choice of the gauge sector is not phenomenologically viable. In 
particular since there is only one 17(1) factor, the heavy Z' particle inherits couplings to the SM 
fermions that are experimentally excluded. Adding a second 17(1) factor without additional model 
building in the exact Z 2 limit is also problematic, since it leads to the existence of a second mass¬ 
less photon. A number of phenomenological studies of the Twin Higgs model and its variants have 
focused on these and other issues [301(69] however for the purposes of this paper we choose to 
work with this very minimal model. While extended models exist that address such phenomeno¬ 
logical issues, using such a model would only obscure the simple idea behind our analysis without 
significantly altering our conclusions^} 

The cancellations to the Higgs mass arise from an approximate 517(4) global symmetry in the 
scalar sector, of which the SU(2)a x SU(2)b subgroup is gauged. The scalar degrees of freedom 
belong to the fundamental representation of this global SU{ 4) symmetry, such that under the 
gauged subgroup they transform as 


H a 


SU(2) a 


Hr 



SU(2) b 


H a 

H b 


( 2 ) 


Up to a term that will be added later, the tree level potential for the scalars is chosen to respect 
the global 517(4) symmetry, 



(3) 


The SU (4) symmetry is spontaneously broken down to SU (3) as H acquires a vacuum expectation 
value (VEV), which results in seven Nambu-Goldstone bosons and a heavy radial mode. Below the 
scale /, the radial mode can be integrated out to obtain a nonlinear sigma model for the degrees 


Most extended models need to introduce additional breaking of the Z 2 symmetry, and deviations from the exact 
symmetry limit tend to reintroduce quadratic divergences which lead to O ( T 2 ) symmetry restoring mass terms. 
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FIG. 1. The cancellation of quadratic divergences in the top sector in terms of Feynman diagrams in the 
linear (top row) and non-linear (bottom row) formulations of the model. 


of freedom parameterized as 


exp 


/ 



hi 

0 

h 2 


h3 

h\ h\ hi 

ho 


if 




(4) 


which defines h as the SM Higgs doublet field, and h! as the twin Higgs which is charged under 
the twin SU( 2). It is straightforward to see that to leading order, h = Ha- 

The global SU(4) is broken down to SU(2) A x SU(2) B when the theory is gauged, and once 
the Yukawa interactions are introduced, where Ha couples Qa and T A , and Hb couples Qb and 


L B- 


C Yukawa = V ( H\Q A T% + H%Q B T%) 


(5) 


Note that these terms are compatible with the Z 2 symmetry even though they explicitly break 
the SU(A). This has a very important consequence: one-loop corrections to the quadratic part of 
the scalar potential respect the Z 2 , that is, they are proportional to H^ a Ha + H B H B , which can 
be written as In other words, leading quantum corrections to the quadratic part of the 

potential accidentally respect the full global SU{A) symmetry. Specifically, corrections from the 
Yukawas and the SU (2) gauge groups have the following form: 


Hi (H) D 


3y 2 A 2 9g 2 A 2 


h\h a + H' b H b 


( 6 ) 


87t 2 647t 2 

Therefore, any quadratically divergent contributions give a mass to the radial mode of H^H, but 


In fact, quadratically divergent mass corrections have this property to all loop orders. Even if the Z 2 symmetry 
is softly broken by the term to be introduced later in this section, there will be higher-loop mass corrections 
proportional to /r 2 but those are not quadratically divergent. 
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not the SM Higgs doublet h. This cancellation is easiest to see from the linear theory, and appears 
to be somewhat mysterious from the point of view of the low energy theory due to an unusual 
four-point coupling between the SM Higgs doublet and the partner fermions. This is illustrated 
in figure [lj The Z 2 similarly prevents divergent contributions from the gauge sector, which again 
is most easily seen in the linear theory, but of course this also holds true in the nonlinear sigma 
model of the low energy theory after the radial mode has been integrated out. 

More quantitatively, the one-loop Coleman-Weinberg potential (in Landau gauge) 

includes contributions from the top sector, where 


Vcw[H) = 64^ STr 


2 2 f 2 ■ 2 ( v \ 1 2 2 f2 2 ( v \ 

m tA =y f sin J and rn tR =y f cos j , 

with (h) = ^ (u,0), and from the gauge sector, where, in the gjj{i) Y 0 limit 


( 8 ) 


2 f2 


m W A = 


9 Z f 


sin 


V2 f 


_ 9 2y lw j 2 9 2 f 2 

= — and m WB = cos 


_,2 f 2 


V2 f 


(9) 


Since the tree level potential thus far respects the SU (4) symmetry, no potential for h is generated 
from the scalar sector at one loop. Thus h only acquires a mass at one-loop through the top and 
gauge sectors, with the former dominating over the latter. 

The scale of electroweak symmetry breaking uew is defined in terms of the gauge boson masses, 
as shown in eq. <©• Since the Higgs particle is among the non-linearly parameterized Goldstone 
modes, y/2 (h) = v ^ vew = 246 GeV (see figure^. As discussed in ref. [59], this implies that the 
coupling of the Higgs to the weak bosons would deviate from the SM predicted values by a factor of 
cos ■ Exact Z 2 symmetry implies uew = / and that the Higgs couples with equal strength to 

both A and B sector gauge bosons. For this reason, exact Z 2 symmetry is not phenomenologically 
viable. 

If we assume that the exact Z 2 is broken such that wew ^ f ■> with the partner sector being 
heavier than the SM, the mass of the Higgs is set roughly as 


2 3 y 2 2 

m h ~ ^2 m t B lQ g 


A 2 


( 10 ) 


So, for mh = 125 GeV, we are led to expect / ~ 500 GeV, which also justifies the assumption 
of Z 2 breaking. For more details on phenomenological considerations in Twin Higgs models and 
experimental consequences, see ref. [301469] 














7 


H a 



FIG. 2. Graphical illustration of parameters in field space: / and v ^ «ew■ 

To achieve a soft breaking of the Z 2 symmetry, a term /j, 2 H a Ha is added to the potential. Note 
that having a small /x 2 is technically natural since it is the only coupling in the theory that violates 
the Z 2 symmetry. So, for /x 2 ~ ?n 2 (which itself arises at one-loop), higher-loop effects of /i 2 can 
be safely neglected. 


III. THE EFFECTIVE POTENTIAL AT FINITE TEMPERATURE 


In this section we review the basic aspects of finite temperature field theory, which we need 
to compute the effective potential for the Twin Higgs model in the finite temperature equilibrium 
state. We use the Matsubara formalism for finite temperature calculations EDI HQ. 

Let us consider a renormalizable field theory in the perturbative regime, where we only turn 
on a background value for one scalar degree of freedom, denoted from here on as 4>. Since we are 
interested in the phase structure of a gauge theory in particular, the scalar in question will be 
taken to transform nontrivially under a gauge group. The full one-loop finite temperature effective 
potential for (f> can be split up into a zero temperature part (including one loop effects) and a finite 
temperature correction 

14ff(</>, T) = U tree (0) + Lf =o (0) + AV? (</>, T) (11) 


with 


ALf(^,r) 





( 12 ) 


where for each particle denoted by the label i , mj(</>) denotes its mass in the background (j>, and 
by our assumption of perturbativity < O ( 1 j)). The supertrace includes the correct factor 

accounting for the number of degrees of freedom associated with each particle and a minus sign 
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for fermions. Jb and J$ arise from the Bose-Einstein and the Fermi-Dirac distribution functions 

nrr\ 2 ( 


fC 

— / 

Jo 


dt t 2 log 


as 

-y/xi+t 2 

(13a) 

-yj Xi+t 2 

(13b) 


poo 

Jf( x i)= / dtt 2 log 
Jo 

Note that due to the gauge symmetry, any phase of <f> is equivalent, and from this point on we will 
restrict ourselves to 4> > 0. 

While the effective potential is not a gauge invariant object, the value of the potential in the 
vacuum state is well-defined UMM- Since we are only interested in the question of whether the 
symmetry is broken, rather than the details of the phase transition, we can simply investigate 
whether the global minimum of the finite temperature effective potential occurs at the origin of 
field space, defined as the point where the gauge bosons are massles^J 

For any given value of 0, we will mainly be interested in high temperatures T 2 > </> 2 which 
due to perturbativity is equivalent to T 2 > m 2 (</>) as mentioned above, and we will often drop 
the subscript to write x < 1 to denote the high temperature regime. In this limit the formulae 
above can be expanded in an asymptotic series (henceforth referred to as the high-temperature 
expansion) 


7T 4 


7r 2 

- 2 

7T£2 x i 

( X \ 


(14a) 

” 45 

+ 

12 X ~ 

6 32 log 1 


+ ... 

7vr 4 


7T 2 

X 2 1 

^ i ) 


(14b) 

" 360 

+ 

24 ' X 

+ 32 log l 

V®// 

| +... 


where ap = 7r 2 e _27B+ ^ and a& = 167r 2 e~ 27E+ ^. 

In figure [3] we compare, for bosons and fermions respectively, a numerical evaluation of equa¬ 
tions (fl3|) to the truncation of equations (14) at linear order for x, and to a truncation up to 


and including the logarithmic terms. Inspecting the figure, it is evident that the O (x) truncation 
captures the one-loop effective potential only at very high temperatures (T > m), while the log(x) 
truncation does so at roughly T > m or even slightly lower temperatures. 


Let us consider the salient features of the high-temperature expansion in equation (14) term by 
term, starting with the largest thermal contributions. 

Terms of O : Both bosons and fermions have ^-independent Stefan-Boltzmann contribu¬ 
tions ~ T 4 . This does not affect the structure of symmetry breaking. 


To be precise, the point where the transverse polarizations of the gauge bosons are massless, at the perturbative 
level. 


4 
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T/m T/m 


FIG. 3. Comparison of different truncations of the high temperature effective potential: fermions to the left 


and bosons to the right. Solid (red) lines represent the numerical evaluation of eq. (13), dotted (green) lines 
and dashed (blue) lines respectively represent truncations to linear order in x and a truncation up to and 


including the logarithmic terms, in eq. (14). 


Terms of O (x 1 ): This is the first order at which V e R picks up a (f> dependence, and since 
T 4 x = T 2 m 2 ((j)) ~ g 2 T 2 (j) 2 where g symbolically denotes the strength of coupling between 0 and 
the particle labeled by the index i, the contributions of O (x) provide an effective thermal mass 
for </> proportional to gT (masses in the EFT are parametrically smaller than the Matsubara scale 
7 rT). At finite temperature, bosonic and fermionic modes running in a loop contribute to this term 
with the same sign because they have opposite boundary conditions on the thermal circle. This 
is connected to the fact that supersymmetry is broken at finite temperature, and the scalar mass 
term can acquire large positive corrections <5m 2 h ~ g 2 T 2 in a supersymmetric theory even though 
contributions to m 2 cancel at zero temperature. These contributions to the effective thermal mass 
of <j) generically drive the scalar background value towards the origin of field space. For bosons and 
fermions respectively, one can set up a correspondence between thermal mass corrections and zero 
temperature divergent mass contributions m 



A 2 

T 2 


bosons: 

167T 2 

12 

(15a) 


A 2 

rp2 


fermions: 

167T 2 

24' 

(15b) 


On the other hand, a symmetry that leads to cancellations between the contributions of same- 
spin particles to the mass of ^ at zero temperature will also induce a cancellation among corre¬ 
sponding thermal mass contributions, which makes symmetry non-restoration at finite temperature 
a possibility. This is precisely the case in models where the Higgs is embedded as a pNGB in a 
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FIG. 4. Schematic diagram of daisy resummation. Blue (dashed) lines correspond to zero-modes and the 
blobs correspond to the loops involving higher Matsubara modes. 

nonlinear sigma model and the coupling of the Higgs to the heavy fermionic partners arises from 
higher dimensional terms. This is how the divergences can conspire to cancel at zero temperature, 
as illustrated in figure [1} This cancellation can be preserved when the model is UV-completed into 
a linear sigma model, which is true in the linear formulation of the Twin Higgs model presented 
in section El 

If the O (x) terms can be made to cancel in this fashion, then the phase structure of the model 


will depend on the effect of the subleading terms in equation (14) which therefore must not be 


neglected. The physics behind these terms is more subtle and we discuss them next. 

Terms of O (x 3 / 2 ): In the Matsubara formalism one can expand the fields into their Kaluza- 
Klein modes around the compact thermal direction. All heavy modes can then be integrated out, 
leaving us with a dimensionally reduced effective field theory (EFT) of the zero modes, in three 
(spatial) dimensions. Note that due to their boundary conditions, fermionic degrees of freedom do 
not have zero modes and therefore the EFT is a theory of scalars and gauge bosons only. As can 
be seen in temporal gauge, the gauge boson degrees of freedom arrange themselves into an adjoint 
scalar {A T ) and a gauge field Ai. By dimensional analysis in this EFT, corrections to the vacuum 
energy from zero modes running in loops must be proportional to m 3 (</>), which is nothing but the 


x 3 / 2 term in equation (14). 


Massless zero modes running in loops lead to infrared divergences at higher loop order, which 
means that a one-loop calculation is inadequate close to the phase transition. A better way to deal 
with the zero modes is to include the O (x) thermal corrections to their masses discussed above 
and use the full thermal mass m 2 h (0, T ) = m 2 ((j)) + 5m 2 h (T) when calculating their contribution to 
the one-loop effective potential^} This is equivalent to resumming a series of higher-loop diagrams 
known as “ring diagrams” or “daisies” (illustrated in figure [4]) that capture the most egregious 
infrared divergences, and it is particularly important at temperatures where the thermal mass 


5 As shown in ref. [76] ( A T ) acquires a positive mass, which allows us to restrict our attention to only the 4d scalars 
as the order parameter for the phase transition. 
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correction is comparable to, or larger than m 2 (</>). The ring-corrected finite temperature effective 
potential thus becomes 


Ves(4 h T) = V tiee ^) + Vf=°(^) + AVf (& T) + AF ring (0, T), 


(16) 


with 


Airing = Y, 


12tt 


m 2 {4>) 

rp'2 


m 2 ((j)) + 5m 2 h (T) 

r-p2 


(17) 


where the summation runs over all the scalar degrees of freedom in the dimensionally reduced EFT. 

At high enough temperatures m 2 h ((/>, T) becomes positive, even for scalars which have rri 2 ((f>) < 0 
at zero temperature. This eliminates contributions to the effective potential coming from the x 3//2 
term that naively appear to be imaginary m Note that in the Twin Higgs model all scalar modes 
are pNGBs of the <S[/(4) symmetry at tree level and therefore do not contribute to the one-loop 
potential for the SM Higgs which is also among the pNGBs, so this particular issue does not arise. 

For zero modes of the transverse polarizations of gauge bosons, residual gauge symmetry in 
the dimensionally reduced EFT prevents any perturbative mass corrections, including ring dia¬ 
grams. However, the gauge coupling in the 3d EFT has dimensions of mass, which leads to a 
non-perturbative mass correction Sm 2 p ~ g 4 T 2 [[75]. Lattice results [79] indicate that such non- 
perturbative corrections cannot be neglected in the case of the SM, since they affect the nature of 
the phase transition and reveal the correct expansion parameter of any perturbative description to 

be -7J-. 

m w 

Even using the best analytical methods available, studying the phase transition is a hard prob¬ 
lem [80]. This should not be surprising, since a phase transition corresponds to a non-analyticity 
in how the free energy depends on the parameters of the model, which cannot be captured at any 
finite order of perturbation theory m • While lattice methods are the most reliable approach in 
cases such as second order phase transitions, for stronger phase transitions we can gain a quali¬ 
tative understanding by using analytical methods [.81]. We pursue the latter approach here, and 
hope that our conclusions may serve to motivate further analysis by others. 

Taking into account the corrected masses of gauge bosons in the effective theory, their contri¬ 
bution to the effective potential at high temperature has the following form: 

3 


E T[m 2 ^) + 5m 2 th (T)}i * E ^ 


polarizations 


polarizations 


3 m 2 ((j)) 

1 + + ... 


2 C T 2 


(18) 


E Q^T 4 +-yfQT 2 m 2 ^) + 0{T Q ) 

polarizations 
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where 5rnf h (T) = <(T 2 and ( contains numerical factors and couplings. The T 4 and T 2 terms 
imply corrections to the Stefan-Boltzmann term and the thermal mass of <f> respectively, followed 
by corrections with non-positive powers of T. Strictly speaking, at temperatures where non- 
perturbative thermal mass corrections to the mass of the transverse polarizations of the gauge 
bosons dominate mass contributions coming from the Higgs mechanism, the Higgs VEV ceases to 
be a good order parameter. 

If a cancellation among the O (T 2 ) terms in the one-loop effective potential persists in the three 
dimensional EFT after the resummation, then keeping the subleading terms of O (log T) becomes 
crucial to any attempt at an (approximate) analytical study of the phase structure of the theory. 
For the same reasons as in the discussion of the O (x) terms, this is indeed the case for the model 
at hand so we finally turn our attention to this last set of terms. 

Terms of O(logx): Note that for each degree of freedom, the logarithmic terms in eq. © 
combine with the logarithmic terms in the zero temperature Coleman-Weinberg potential of eq. (Jt[) 
to give a log-^ 5 — dependence on the temperature, as the factors of m 2 ((j)) cancel between the 
one-loop corrections at zero and finite temperature. Any formal cutoff dependence thus comes 
from the zero temperature Coleman-Weinberg potential, whose parameters have been chosen to 
reproduce the observed electroweak VEV and Higgs mass. As we will see in the next section, 
the non-cancellation of these terms will determine the fate of electroweak symmetry restoration at 
finite temperature, in the Twin Higgs model we consider. 

Having reviewed the most important aspects of field theory at finite temperature in general, we 
will apply what we have learned specifically to the Twin Higgs model in the next section. 


IV. TWIN HIGGS AT FINITE TEMPERATURE 

Let us now specialize our discussion to the Twin Higgs model at finite temperature, and let 
us consider whether there can be any important contributions to the effective potential that we 
have not already accounted for in the previous section. Due to invariance under gauge symmetries, 
H^Ha and H^Hb are the only combinations that the effective potential can depend on at zero 
or finite temperature, and in the exact Z 2 limit these have identical coefficients, which ensures 
an accidental S'C/(4) symmetry for quadratic terms, forbidding any dependence on the Goldstone 
modes at that order. Furthermore, even though the Z 2 is broken by the /r 2 term introduced at 
the end of section [Hj this is a soft breaking, thus any Z 2 violating corrections to the potential 
must include a positive power of /r 2 . This means that by simple dimensional analysis, at the 
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FIG. 5. We plot the finite temperature effective potential at two different temperatures: The blue solid line 
represents the potential (numerically evaluated) at Ti = 100 GeV, the red dashed line represents the potential 
(numerically evaluated) at T 2 = 350 GeV and the green dotted line represents the O (T 2 ) truncation of the 
potential. Note that for the quadratic truncation, the potential is independent of temperature, and hence, 
does not sense symmetry restoration. See the main text for the numerical values of the relevant parameters 
that were used in making the plot. 


renormalizable level there can be no Z 2 violating contributions to the potential with a positive 
power of temperature that depend on the Goldstone modes either. There can be contributions of 
order log(T) that are Z 2 violating, but these are subdominant to the contributions of order log(T) 


that have already been considered at the end of section III, since «/"• 


Of course, the potential has a dependence on the Goldstone modes beyond quadratic order, 
where Z 2 invariance is no longer equivalent to full SU (4) invariance, leading to terms such as 
+ (^H^Hb^■ Being dimension four operators, by dimensional analysis the coefficients 
of such terms also cannot include a positive power of temperature, so these fall into the class of 
contributions of order log(T) that we have already discussed. 

Since we have now convinced ourselves that all thermal mass corrections of O (T 2 ) cancel, let 








14 


us proceed to evaluate the logarithmic contributions from the top quark and its partner: 

2 / , \ 2" 


- (3 x 4) 


{y 2 f 2 sin 2 j) (y 2 f 2 cos : 


2 h 

f 


647r 2 


647T 2 


log 


a F T 2 

A 2 


= -12 


(y 2 {/ l + ...} 2 ) 2 + (y 2 { /_ ^ + ---} 2 ) 


647T 2 


647T 2 


log 


cl f T 2 


(19) 


3y 4 f 2 h^h a F T 2 

In hindsight, the fact that these subleading contributions do not cancel each other should come 
as no surprise. The non-cancellation of logarithmic terms between Z 2 partners in the one-loop 
Coleman-Weinberg potential is precisely what keeps the Higgs from being an exact Goldstone 
boson and ensures a sizable Higgs mass at zero temperature (see for instance equation (|10[) ). As 
mentioned previously, note that the appearance of A eq. [19] arises from the zero temperature 
Coleman-Weinberg potential, and the finite temperature additive corrections are independent of 


the cutoff (refer eq. 14). 

More generally, in the case of phenomenologically viable SM extensions with a pNGB based 
mechanism for naturalness, any ultraviolet divergent contribution to the Higgs effective potential 
at zero temperature will carry over to a corresponding finite temperature contribution. In the 
specific case of the Twin Higgs model, as well as in other natural models with a similar symmetry 
structure, we expect this feature to drive symmetry restoration at finite temperature. 

If we wish to study the phase of the theory at temperatures around the electroweak scale (and 
not significantly higher than the partner masses), it is straightforward to numerically evaluate the 


one-loop effective potential. While the high-temperature approximations in eq. (14) are analytically 


tractable and help shape our thinking, we choose to numerically evaluate eq. (12) and eq. (13) in 


order to avoid any artifacts from truncating the expansion. The results for our benchmark model 
(with / = 450 GeV, = 90 GeV and A = 4.4 TeV) are presented in figure [5] The most important 
one-loop effects come from the top sector, followed by the electroweak gauge sector, resulting 
in the restoration of electroweak symmetry at T ~ 300 GeV (which stays restored as we push 
temperatures up to where the EFT starts breaking down). 

Since non-perturbative effects cloud the study of physics close to the phase transition, it would 
be nice to attack this question from a different angle. In particular, it would be of interest to 
look for a symmetry restored phase at temperatures much higher than the phase transition, where 
a resummed theory has a valid perturbative description. Of course, this cannot be done in the 
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nonlinear Twin Higgs model and necessitates a UV-completion, which we take to be a linear 
sigma model completion of the Twin Higgs described in section [TTJ In this UV-completion, the 
“radial mode” linearizes the sigma model. The radial mode is a singlet of the approximate SU (4) 
global symmetry in the scalar sector, and thus its zero-temperature mass will not be protected from 
quadratically divergent contributions. This means that at high temperatures, the radial mode picks 
up a positive thermal mass term, leading to a symmetry restored phase in the UV completion. The 
radial mode being driven to zero is a sufficient condition for the gauge bosons to become massless 
(up to thermal contributions). Note that this result is not directly related to our calculation in 
the nonlinear model, as the VEVs of the radial mode and the Goldstone modes are separate from 
each other. 


V. CONCLUSIONS AND OUTLOOK 

We have investigated the possible existence of a broken phase of electroweak symmetry at high 
temperature in extensions of the Standard Model where the Higgs is realized as a pNGB, focusing 
on the Twin Higgs model as a benchmark. While we have confirmed that one-loop quadratic 
contributions to the Higgs potential at finite temperature cancel between the Standard Model 
degrees of freedom and their partners as they do at zero temperature, this is not true for subleading 
corrections to the effective potential, which restore electroweak symmetry at high temperature. 
Cancellation of O (T 2 ) corrections to terms in the the Higgs potential is a generic consequence 
of same-spin partners ensuring naturalness at zero temperature, and the logarithmic corrections 
are connected to obtaining a phenomenologically viable Higgs boson mass at zero temperature. In 
the case of the Littlest Higgs model considered in ref. [7], the EFT has uncancelled quadratically 
divergent corrections to higher-order terms in the Higgs potential (arising from non-renormalizable 
operators), but even in that case, the theory exhibits a restoration of electroweak symmetry as long 
as temperatures are not pushed beyond the range of validity of the EFT for a finite temperature 
calculation. 

It should be noted that nonlinear sigma models in which the Higgs is a pseudo-Goldstone 
generically become strongly coupled at high energies and require UV completion, which brings up 
the question of whether a suitable UV completion may nevertheless allow for a broken phase of 
electroweak symmetry to persist at high temperature. As we have demonstrated for the case of the 
Twin Higgs model, UV completing the theory into a linear sigma model cannot achieve this, since 
the mass of the radial mode is unprotected from quadratic corrections, which at finite temperature 
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drive the radial mode to the origin of field space and lead to symmetry restoration. One can also 
contemplate nesting one nonlinear sigma model inside another with a higher symmetry breaking 
scale, however since the Higgs receives a thermal correction (albeit at the subleading level) in the 
original nonlinear sigma model, this type of construction will not change the finite temperature 
behavior. 

Alternatively, one can imagine supersymmetrizing the linear sigma model, since supersymmetry 
is the best understood UV-complete mechanism to protect the mass of a scalar from quadratically 
divergent corrections. However we know that supersymmetry does not prevent quadratic mass 
corrections at finite temperature, and therefore the radial mode VEV would still be driven to zero. 

Let us also briefly remark on classes of natural extensions of the SM other than supersymmetry 
and Higgs as a pNGB. Theories with strongly coupled Higgs sectors appear to be disfavored in 
light of the experimental findings at the LHC, and in any case these typically exhibit symmetry 
restoration for temperatures above the formation of the condensate. Gauge-Higgs models have been 
shown to lead to a restored symmetry phase at high temperature [82]. In “Relaxion” models [83] , 
electroweak symmetry is also restored at high temperatures, but this idea is quite recent and it 
would be interesting to study whether variants of it may have a more subtle finite temperature 
behavior. 

Our conclusions are also consistent with more general considerations based on the thermody¬ 
namic behavior of systems at high temperature. In particular, the free energy of a system is given 
by 


T = £ — TS, 


( 20 ) 


and therefore, heuristically, at high temperatures, the free energy can be minimized by increasing 
entropy (corresponding to a symmetric phase) rather than lowering the energy by spontaneous 
symmetry breakdown [83]. This suggests a robust rule of thumb that symmetries get restored at 
high temperatures, in the absence of any other thermodynamic variables describing the system that 
can attain values that are “natural” based on dimensional analysis. We remark in passing that 
if this last criterion is removed, e.g. when the system has a chemical potential /r ~ T, symmetry 
non-restoration is possible, see for example refs. |85H88| . 
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